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Abstract 

We study the stability of circular orbits of the electromagnetic two-body problem in an electro- 
magnetic setting that includes retarded and advanced interactions. We give a method to derive 
the equations of tangent dynamics about circular orbits up to nonlinear terms and we derive the 
linearized equations explicitly. In particular we study the normal modes of the linearized dynamics 
that have an arbitrarily large imaginary eigenvalue. These large imaginary eigenvalues define fast 
frequencies that introduce a fast {stiff ) timescale into the dynamics. As an application of Dirac's 
electrodynamics of point charges with retarded-only interactions, we study the conditions for the 
two charges to perform a fast gyrating motion of small radius about a circular orbit. The fast 
gyration defines an angular momentum of the order of the orbital angular momentum, a vector 
that rotates in the orbital plane at a frequency of the order of the orbital frequency and causes 
a gyroscopic torque. We explore a consequence of this multiscale solution, i.e; the resonance 
condition that the angular momentum of the stiff spinning should rotate exactly at the orbital 
frequency. The resonant orbits turn out to have angular momenta that are integer multiples of 
Planck's constant to a good approximation. Among the many qualitative agreements with quan- 
tum electrodynamics (QED), the orbital frequency of the resonant orbits are given by a difference 
of two eigenvalues of a linear operator and the emission lines of QED agree with our predictions 
within a few percent. 
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I. INTRODUCTION 



We study the Lyapunov stability of quasi-circular orbits of the electromagnetic two-body 
problem, a dynamical system with implicitly-defined delay. The motivation is to understand 
the balancing of the fast dynamics described by the delay equations for particle separations 
in the atomic magnitude. We give an economical method to derive the linearized equations 
of motion about circular orbits. We work in a generalized electromagnetic setting where the 
field of the point charge is a linear combination of the advanced and the retarded Lienard- 
Wiechert fields with an arbitrary coefficient [1], henceforth called the Eliezer setting (ES) 
(see Appendix A) . We study in detail a specific feature of the tangent dynamics; The stiff 
normal modes of the linearized dynamics, that have an arbitrarily large imaginary eigen- 
value and introduce a fast timescale in the dynamics. The derivation of the fast normal 
modes of tangent dynamics is our main technical contribution to the understanding of this 
dynamical system. For some special cases of the ES we find a remarkable quasi- degeneracy 
of the tangent dynamics; these cases are Dirac's theory with retarded-only fields, the action- 
at-a-distance electrodynamics and the dissipative Fokker theory of Ref. [2] (see Appendix 
B). Last, we discuss the dynamics of the hydrogen atom in Dirac's electrodynamics with 
retarded-only fields, the special case of the ES of greatest relevance to physics. Having 
recognized the existence of the fast dynamics near circular orbits, our method to find the 
trajectory starts by balancing the fast dynamics near a tentative circular orbit. We inves- 
tigate the conditions for the dynamics to execute a fast gyrating motion about a circular 
orbit, henceforth called spinning and illustrated in Fig. 1. This fast gyration defines an an- 
gular momentum vector of the order of the orbital angular momentum of the unperturbed 
circular orbit. This angular momentum of spinning is a vector that rotates in space with 
a frequency of the order of the orbital frequency. We experiment with a typical necessary 
condition for such multiscale solution; the resonance condition that the angular momentum 
of spinning rotates exactly at the orbital frequency. This resonance condition turns out to 
be satisfied precisely in the atomic magnitude. We predict several features of the hydrogen 
atom of quantum electrodynamics (QED) [3] with good precision and qualitative detail. The 
frequencies of the resonant orbits agree with the lines of QED within a few percent average 
deviation. There is also a large body of qualitative agreement with QED; (i) the emitted 
frequency is given by a difference of two linear eigenvalues (the Rydberg-Ritz principle) (ii) 
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the resonant orbits have angular momenta that are approximate multiples of a basic angu- 
lar momentum. This basic angular momentum agrees well with Planck's constant and (iii) 
the angular momentum of the fast spinning dynamics is of the order of the orbital angular 
momentum. 

Dirac's 1938 work [1] on the electrodynamics of point charges gave complex delay equa- 
tions that were seldom studied. Eliezer generalized Dirac's method of covariant subtraction 
of infinities in 1947[4] to include advanced interactions naturally in the electrodynamics of 
point charges (the ES) [4]. Among the few dynamics of point charges investigated, another 
early result of Eliezer [5-8] revealed a surprising result (henceforth called Eliezer's theorem); 
An electron moving in the Coulomb field of an infinitely massive proton can never fall into 
the proton by radiating energy. The result was generalized to motions in arbitrary attractive 
potentials [6], as well as to tridimensional motions with self-interaction in a Coulomb field 
[7, 8], finding that only scattering states are possible. Since our model has the dynamics 
of Eliezer's theorem as the infinite-mass limit, a finite mass for the proton is essential for a 
physically meaningful dynamics; If the proton has a finite mass, there is no inertial frame 
where it rests at all times, and this in turn causes delay because of the finite speed of light. 
A finite mass for the proton is what brings delay into the electromagnetic two-body dynam- 
ics, with its associated fast dynamics. This infinite-mass limit is a singular limit, because 
the equations of motion pass from delay equations to ordinary differential equations! Our 
understanding of this two-body dynamics might prove useful for atomic physics and perhaps 
we can understand QED as the effective theory of this complex stiff dynamics with delay. 
We shall describe the two-body motion in terms of the familiar center- of-mass coordinates 
and coordinates of relative separation, defined as the familiar coordinate-transformation that 
maps the two-body Kepler problem onto the one-body problem with a reduced mass plus a 
free-moving center of mass. We stress that in the present relativistic motion the Cartesian 
center-of-mass vector is not ignorable, and it represents three extra coupled degrees-of- 
freedom. We introduce the concept of resonant dissipation to exploit this coupling and 
the many solutions that a delay equation can have. Resonant dissipation is the condition 
that both particles decelerate together, i.e., the center-of-mass vector decelerates, while the 
coordinates of relative separation perform an almost-circular orbit, despite of the energy 
losses of the metastable center-of-mass dynamics. Last, we stress that our point charges 
are not spinning about themselves, but rather gyrating about a center that is moving with 



3 



the circular orbit, as illustrated in Fig. 1, even though we refer to this gyrating motion as 
spinning. 

The road map for this paper is as follows; In Section IV we give the main technical part of 
the paper; We outline our economical method to derive the tangent dynamics of the circular 
orbit by expanding the implicit light-cone condition and the action to quadratic order. This 
economical method shall be useful in the further research needed to derive the higher orders 
and the complete unfolding of this dynamics. In this Section we also take the stiff limit 
of the linear modes of tangent dynamics and introduce the fast timescale. In Section V 
we give an application to atomic physics, by discussing a consequence of balancing the fast 
dynamics first; This balancing defines an angular momentum of fast gyration comparable 
to the orbital angular momentum, a vector that rotates in space with a slow frequency of 
the order of the orbital frequency. In this section we discuss the heuristic condition that 
the angular momentum of the fast spinning motion should rotate at the orbital frequency, a 
resonance condition that predicts the correct atomic scales. The earlier sections are a prelude 
to Section IV. In Appendix A we give the equations of motion of the electrodynamics of point 
charges in the ES. Section II is a review of the circular orbit solution, to be used in Sections 
IV and V. In Section 111 wc give an action formalism for the Lienard-Wiechert sector of the 
ES as a prelude to the quadratic expansions needed for the linear stability analysis of Section 
IV. In Appendix B we derive the general tangent dynamics for oscillations perpendicular to 
the orbital plane. The Lemma of resonant dissipation of Appendix C proves that for the 
linearized equations of motion the state of resonant dissipation is impossible. In Appendix 
D we discuss how the nonhnear stiff terms balance the leading dissipative term of the self- 
interaction force and we estimate the radius of the stiff torus. Last, in Section VI we put 
the conclusions and discussion. 

II. THE CIRCULAR ORBIT SOLUTION 

The circular-orbit solution of the isolated electromagnetic two-body problem of the action- 
at-a-distance electrodynamics (the ES with k = —1/2, see Appendix A) is used here as the 
unperturbed orbit[9, 10]. For the ES with k — —1/2, the tangent dynamics of the next 
section is straightforward Lyapunov stability analysis. In the general case, ( A; 7^ —1/2), the 
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ES also prescribes a small force opposite to the velocity of each particle along the circular 
orbit, such that the circular orbit is not an exact solution of the equations of motion. The 
tangent dynamics is then the starting point of a perturbation scheme to obtain a stiff torus 
near a circular orbit (the state of resonant dissipation discussed in Appendix C). 

We shall use the index i = 1 for the electron and i = 2 for the proton, with masses mi and 
777.2 respectively. We henceforth use units where the speed of light is c = 1 and ei = —62 = — 1 
(the electronic charge). The circular orbit is illustrated in Fig. 2; the two particles move in 
concentric circles with the same constant angular speed and along a diameter. The details 
of this relativistic orbit will be given now; The constant angular velocity is indicated by 
Q , the distance between the particles in light-cone is and 6 = Qr^ is the angle that one 
particle turns while the light emanating from the other particle reaches it (the light-cone 
time lag). The angle 9 is the natural independent parameter of this relativistic problem. For 
orbits of the atomic magnitude the orbital frequency is given to leading order by Kepler's 
law 

n^iie^ + ... (1) 

while the light-cone distance rb is given by 

r. = ^ + ... (2) 

Each particle travels a circular orbit with radius and scalar velocity defined by 

ri = bin, (3) 
r2 = 62^6, 

vi = = ebi, (4) 

for the electron and for the proton, respectively. The condition that the other particle turns 
an angle 9 during the light-cone time lag [9] is 

bl + bl + 26162 cos(^) = 1, (5) 

and is henceforth called the unperturbed light-cone condition. As shown in Appendix B 
of Ref. [2], bi and 62 are approximated by the Keplerian values 61 = 1112/ {mi + 777.2) and 
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^2 = mi/ (mi + 1712) plus a correction of order 9^. As discussed in Ref. [9], there is a 
conserved angular momentum perpendicular to the orbital plane of the circular orbit given 

by 

_ 1 + t;it;2COs(6') 
'~ 9 + v,V2sm{9)' ^ ^ 

where the units of are e^/c, just that we are using a unit system where — c — 1. For 
small values of 9 the angular momentum of Eq. (6) is of the order of Iz ~ 9~^. For orbits of 
the atomic magnitude, Iz — 9~^ is about one over the fine-structure constant, — 137.036. 



III. AN ACTION FOR THE LORENTZ FORCE 



We introduce an action to be used as an economical means to derive the Lorentz-force 
sector of the ES equations of motion with the linear combination of retarded and advanced 
Lienard-Wiechert potentials. In the following we derive the tangent dynamics by expanding 
this action to quadratic order. We henceforth use the dot to indicate the scalar product of 
two Cartesian vectors. The Lienard-Wiechert action for the Lorentz-force sector of the ES 
is 

e = -kf + + ('-v.- v..) 

J ri2a(l-Fni2a • V2a) J ri2b{l - n^b ' V2b) 

In Eq. (7) Vi stands for the Cartesian velocity of particle 1 at time ti and V2a and V26 
stand for the Cartesian velocities of particle 2 at the advanced and at the retarded time t2 
respectively. The vector ni2a is a unit vector connecting the advanced position of particle 2 
at time t2 to the position of particle 1 at time ii, vector ni26 is a unit vector connecting the 
retarded position of particle 2 at time t2 to the position of particle 1 at time ti and (Below 
Eq. (8) of Ref. [2] this normal is defined to the opposite direction, but it was used correctly 
in Eq. (19) of Ref. [2] ). Notice that each integral on the right-hand-side of Eq. (7) can be 
cast in the familiar form 

' A^-^VI^^^i ^- f{V-Vi- A)dt, (8) 

where V and A are the Lienard-Wiechert scalar potential and the Lienard-Wiechert vector 
potential respectively. We have introduced the quantity c = ±1 in the denominator of Eq. 
(8) such that c—1 represents the advanced interaction while c——l represents the retarded 
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interaction. The quantities of particle 2 in Eq. (8) are to be evaluated at the time t2 defined 
implicitly by 

t2^h + 'f, (9) 

where c—±l describes the advanced and the retarded light cones, respectively. Because of 
this decomposition into V and A parts, we henceforth call Eq. (8) the VA interaction. The 
minimization of action (7) plus the kinetic energy of particle 1 yields the equations of motion 
for particle 1 suffering the Lorentz-force of the Lienard-Wiechert potentials produced by the 
other particle (see Eq. (58) ), as shown for example in Ref.[ll]). 

Besides the Lorentz-force, which is the Lagrangian part of the ES equations of motion, 
there is also the dissipative self-interaction force in Eq. (58). The shortest way to write the 
ES equations of tangent dynamics is to add this term to the Lagrangian equations of motion, 
watching carefully for the correct multiplicative factor. The stiff limit is determined by the 
largest-order derivative appearing in the linearized equations of motion of Appendix A. In 
this approximation, the contribution of the self-interaction force of the ES to the linearized 
dynamics about a circular orbit is simply given by the Abraham-Lorentz -Dirac force 

Fracl = l{l + 2k)k, (10) 

with a renormalized charge. The contribution of the other smaller nonlinear terms will be 
given elsewhere (one such nonlinear term is used in the estimate of Appendix D). 



IV. EXPANDING THE ACTION TO FIRST ORDER: MINIMIZATION 



In the following sections we substitute the circular orbit of Eqs. (3) and (4) plus a planar 
perturbation into Eliezer's action (7). We then expand the action up to the quadratic order 
to yield the linearized equations of motion. Our economical derivation of the stability 
equations anticipates the need to study the stability of the stiff torus near the circular orbit, 
a novel solution of the two-body problem as illustrated in Fig 2. In this work we keep to the 
linear stability of circular orbits, which is done with the quadratic expansion of the action. 
The variational equations for planar perturbations arc decoupled from the equation for 
transverse perturbations. We perform the planar stability analysis using complex gyroscopic 
coordinates rotating at the frequency fl of the unperturbed circular orbit; The coordinates 
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{xj,yj) of each particle are defined by two complex numbers rjj and { j — 1 for electron 
and j — 2 for proton) according to 

Uj = Xj + iyj = rbexp{iQt)[dj + 2r]j\, (11) 
u* = Xj - iyj = r6exp(-int)[(i* + 2^^], 

where di = bi and ^2 = — &2 are defined in Eqs. (3). Because Xj and yj are real, we must 
have rjj = ^* but to obtain the variational equations it suffices to treat rjj and as two 
independent variables in the Lagrangian. We henceforth call the coordinates of Eq. (11) the 
gyroscopic f] and ^ coordinates. The rf's are real numbers if we choose the origin of times at 
t — 0, but we have kept the harmless star above them to indicate complex conjugation in 
anticipation to what follows. Two quantities appear so often in the calculations that we have 
named them; (i) The numerator of the action (8) along circular orbits, henceforth called C, 
evaluates to 

C = 1 + 6162^^ cos(^), (12) 

for either retarded {c — —1) or advanced interactions { c— 1), and (ii) The denominator of 
action (8) along circular orbits, divided by r;, , henceforth called S and defined by 

5 = 1 + 6iMsin(^), (13) 

again for either retarded or advanced interactions {c — ±1). For the stiff limit of Sections 
V and VI we shall ignore the 0{9^) corrections and set C = 5" = 1. Here we shall derive 
the electron's equation of motion only (j = 1 in Eq. (11)). The equation for the proton 

is completely symmetric and can be obtained by interchanging the indices. One can derive 
from Eq. (11) that the electron's velocity at time ti is 

ill = + iviy = 9 exp{iflti)[idi + 2i{r]i — if]i], (14) 
ill = vix — iviy = 9 exp{—iD,ti)[—idl — 2i{^i + i^i)]. 

The velocity of the proton at its time ^2 can be obtained by simple interchange of the indices 
1 and 2, remembering that di — bi is defined positive while ^2 = — &2 is defined negative, 
because at the same time the particles are in diametrically opposed positions on the exact 
circular orbit, such that the exchange operation on the d's is di <(=^ — ^2 (see Fig. 2). 
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The coordinates of particle 2 entering in the VA interaction of Eq. (8) are evaluated at 
a time t2 in light-cone with the the present of particle 1. Because the implicit light-cone 
condition has to be expanded and solved by iteration, it is convenient to define a function 

as 

t2^ti + - + |f. (15) 
c iz 

The above definition is good for both the advanced and the retarded cases, c—1 defining 
the advanced light-cone and c = — 1 defining the retarded light-cone and the underscore c 
is to indicate that ip^ is a function of c. If the perturbation is zero then = and we are 
along the original circular orbit, where the light-cone lag is the constant r^. We henceforth 
call t the present time of particle 1 (the electron) and we measure the evolution in terms of 
the scaled-time parameter r = Qi. The implicit definition of ip^ by the light-cone condition 
involves the position of particle 2 at the advanced and the retarded time ^2 as defined by 
Eq. (11), 

U2{t + + = nex.^{iVtt2){d2 + 2ri^{T + + ^J], (16) 
uI{t + ce + (pj = rbexp{-int2)[d*2 + 2^2{^ + c0 + cpj]. 

as well as the velocity of particle 2 at the advanced/retarded position 

U2{t + c9 + (fj = i9 exp{iilt2)[d2 + {^2 - ^^2) + 2y?c(^2 - ^2)], (17) 
u*{t + c9 + p>^) = -i^exp(-zOi2)K + (e2 + 4) + 2<^c(e2 + e2)]- 

The linear stability analysis involves expanding the equation of motion to first order in rjj^ 
and ^j^. , which in turn is determined by the quadratic expansion of the action (7) in r]j^ and 
^f.. We must therefore carry all expansions up to the second order in the 77^ coordinates. For 
example the position vector of particle 2 can be determined to second order by expanding 
the arguments of 77 and ^ in a Taylor series about the value on the unperturbed circular 
hght-cone for one order only as 

U2{t + C0 + ip^) ~ rbexp{iQt2){d2 + 2[r/2(T + c9) + ^pj]2{'T + cO)]), (18) 

uI{t + ce + if,) ~ nexp{-int2){d; + 2[^^{t + cO) + p^i^{T + cO)]}. 
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We shall henceforth indicate the quantities of particle 2 evaluated exactly on the hght-cone by 
writing that quantity with a subindex c, as for example 772c = '72(''' + c^) ^2c = C2(''' + c^)- 
In the following we find that ip^ is linear in 77 and ^ to leading order, such that the next 
term in the above expansion involves third order terms, which are not needed for the linear 
stability analysis of the circular orbit. We shall always expand quantities evaluated at the 
perturbed light-cone in a Taylor series about the unperturbed light-cone up to the order 
needed, such that our method yields equations with a fixed delay. The perturbed light-cone 
is expressed implicitly by the distance from the advanced/retarded position of particle 2 to 
the present position of particle 1, described in gyroscopic coordinates by the modulus of the 
complex number 

Uc = uiir) - U2{t + c9 + (p^). (19) 

where c — 1 describes the advanced position and c — —1 describes the retarded position. 
Using Eq. (18) this complex number evaluates to 

Uc = rbex.p{int2){D* + 2[exp(-ic^ - iiPc)Vi(^) - V2c - 'PcV2c]}^ (20) 
ul = rbexp{-int2){D + 2[exp(ic^ + i(Pc)^i{r) - ^2c " <^cC2c]}> 

where we have factored the ftt2 out and defined the complex number 

D = b2 + bi exp{ic0 + icp^) , (21) 

that carries an implicit dependence on (p^ and will later be expanded to second order as well. 
Notice that the star appears in the unusual first line because we are factoring flt2 out in Eq. 
(20). At cp^ — (the unperturbed circular orbit), the complex number D defined by Eq. (21) 
has a unitary modulus. The last term on the right-hand side of Eq. (20) is a quadratic form 
times exp{iQt2) and in the action it only appears multiplied by a counter-rotating constant 
term, i.e., a quadratic form times exp{—iQt2), such that the product is independent of 
t2- Therefore this quadratic term can be integrated by parts and the Gauge term can be 
disregarded, as in any Lagrangian. We henceforth call this a quadratic Gauge simplification 
and we shall use it in several expressions. This simplification applied to integrate the terms 
in 772c ^2c of Eq. (20) by parts yields 

u = r,eMi^h){D* + 2[(1 - i(^J exp(-ice)77i(r) - (1 - <pM}, (22) 
u* = r5exp(-ini2)P + 2[(l + i(^Jexp(ic^)ei(T)-(l-(^je2c]}- 
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where we have also expanded the exponential of (/? in the coefficient of rj^ and up to the 
linear order in (p^, enough to give the correct quadratic action. Last, the velocity of particle 
2 has the following expansion up to a quadratic Gauge 

U2{t + C0 + if^) ~ i^exp(iQi2)M2 + 2(%, - ir/2e)(l - <Pc)i (23) 

«;(T + C^ + 99j ~ -^^ exp(-^l]t2)[rf; + 2(^20 + ^e2c)(l-<^c)] • 

Using the above quantities, the numerator of the Lienard-Wiechert action (8) can be 
calculated as 

1 1 

/i2 = (1 - • V2c) = (1 - 2^1^2c - 2^1^2c)- (24) 

Last, the denominator of the Lienard-Wiechert action (8) can be calculated as 

/t4 = r + ii- — = rb(l + 0) + ^^ + ^^. (25) 
In Eq. (25) we have introduced the scaled delay function by 

= ce4>. (26) 

To relate 6 to the ^7] perturbations we expand the implicit light-cone condition of the 
perturbed orbit up to the quadratic order 

uu* = |x2(t + - + ^^^e) - = c\h - tf = (r, + ^f. (27) 

c \lc 

The hght-cone condition (27) is simpler when expressed in terms of the scaled defined in 
Eq. (26). The expansion of Eq. (27) up to the second order is the following quadratic form 
in 

C(/.2 + 2S(t> = 2[(feiei - ^2^2) + {hiVi - h2ri2)] (28) 
+2[(62ei - hr]2) exp(ic^) + (62/?! - ^1^2) exp(-ic^)] (29) 
- iiV^ exp(ic^) - ^2^1 exp(-ic^)] 

-2v3[62(C2 + %) + 61^2 cxp( -ic6') -h 6ir/2exp (ic6')]. 

The solution to Eq. (28), to ffist order in the ^77 coordinates is given by 
S(t>{i) = [{biCi - + ihVi - ^2%) + (^2^1 - biQ exp(ice) + (62771 - fci^) exp(-ic^)] (30) 
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As an application of the above, we derive the equations of motion for the circular orbits 
of the action-at-a-distance electrodynamics [9] (the ES with k — —1/2 ). The equation of 
motion for can be calculated by expanding the VA interaction of Eq. (7) to linear order 

e = 1 - {[e'{s - i)s + c^]{h + &2cos(e)) + siesin{e) - cos{e))b2} ^'"-^J'\ (3i) 

where the tilde indicates that we scaled © by the factor Qg — C/{rbS) (the value of © 
along the unperturbed circular orbit). Scaling the kinetic energy with the same factor and 
expanding to first order yields 

f. = =^+'-^^^^(r,,+U (32) 

Such that the effective Lagrangian for particle 1 is 

Li = Ti + © (33) 

Since those are linear functions of the Euler-Lagrange equation for is simply |^ = 
or 

m^hn'jie^s^ = [c^ + e^s{s - 1)] (61 + 62 cos(^)) + s{e sm{e) - cos{e))b2. (34) 

This is Eq. 3.2 of Ref. [9], and the equation for rji is the same condition by symmetry (the 
reason why the circular orbit is a solution). The equation of motion for particle 2 can be 
obtained by interchanging the indices 1 and 2 in Eq. (34), representing Eq. 3.3 of Ref. [9]. 
In the next section we carry the expansion of the action to second order, to determine the 
equations of tangent dynamics. 



V. EXPANDING THE ACTION TO SECOND ORDER 



The next term of expansion (32) of the local kinetic energy of particle 1 can be calculated 
with Eq. (14) to be 

Ti ^ To Q—y^ - \'^^ = Q (^1 + '^i) (35) 
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This kinetic energy of particle 1 has the simple quadratic form defined by two coefficients 



= + (36) 

where Mi = {1 + -fl)mi-f^rb9^S/Cand d = {-fj - 1 )mi'j^9'^S/C. We also need the solution 
of Eq. (28) to second order, = 0(i)+ 4>(2) > where 0(1^ is given by Eq. (30) and 0^2) 
calculated by iteration to be 



S(p{2) = (p{ib2Ci - ibir]2 - bi7]2) exp{ic9) + ip{ibi^2 - ihVi - ^1^2) ^wi-icd) (37) 

-<^&2(l2 + '72)-^C'[</'(l)]' 

+2[Ci»7i - Ci%exp(ic^) - ^2^1 exp(-ic^)] 

Next we expand the numerator and the denominator of action (8), Eqs. (24) and (25), up to 
the quadratic order. We use the expansion of the particle separation, Eq. (20), the velocity 
of particle 1 (Eq. (14)), and the expansion of the velocity of particle 2, Eq. (23). The 
quadratic expansion of the action is cumbersome and must be evaluated with a symbolic 
manipulation software (we have used Maple 9). Integrating by parts and disregarding Gauge 

terms, we can bring this quadratic expansion to the following " normal form for particle 1 

?) 

e, = -'-B{n,t^ - ^^Vi) + Un^.m + Inu^, + ^A^n^? (38) 

Yc ■ . Y* 

+ y (C1C2C - C2cCl) + -Y^Vlh2c - V2cVl) + 

Ac • A* • 

-j{ilV2c - V2cil) + -fiVl^l - ^2Vl) + 

+y C1C2C + Y'^iV2c + yCi?72c + Y^iLc 

Notice that in Eq. (38) the coordinates of particle 2 appear evaluated in either the retarded 
or the advanced unperturbed light-cone, which is indicated by the subindex c. The full 
quadratic expansion of the action is composed of the kinetic energy plus the partial actions 
of Eq. (38) evaluated at the retarded and the advanced light-cones and multiphed by the 

13 



respective Eliezer's coefficient, as in Eq. (7) 

Li = Ti - A;e+ + (1 + A;)e_ (39) 

Notice in Eq. (38) the appearance of the angular-momentum-hke binaries inside parenthesis, 
a normahzation achieved by adding a quadratic Gauge term to the Taylor scries. The 
coefficient of each normal-form binary is obtained by a simple Gauge-invariant combination 
of derivatives, for example the magnetic field coefficient B is given by 

In what follows we evaluate the expansion of the action using a symbolic software (Maple 
9). The linearized Euler-Lagrange equation obtained by minimizing action (39) respect to 
the function is a linear functional involving the coordinates > ''7i > C2+ ) V2+ 1 C2- 
r]2_ as well as their first and second derivatives 

hd^l, Vl,^2+,V2+, ^2-,V2-) = (41) 

-[{Nu + e'G)i, + e'Gl,] + [M^{f|, + 2zfj, - rj,) - UuVi - mi] 

-{K^2+ + R*M-) - + y-L-) + + T-'L-) 

-{P+V2+ + P-V2-) - (A+r/2+ + A-%-) + (^+%+ + ^-%-) 

Minimization of action (39) respect to rji, ^2 ^-nd V2 yields three more linear equations, which 
together with Eq. (41) compose a system of four linear delay equations. To include the 
dissipative self-interaction force into Eq. (41) we substitute 9'^ — r^il^ into the definition of 

Ml above Eq. (37), such that the Euler-Lagrange equation of the kinetic energy (36) respect 
to can be written as 

(1 + ^l)rm^i{S/C)r!n'{fii + 2irj, - rj^) (42) 
= rl{S/C) — {mi'yirbUi). 

On the second line of Eq. (42) we recognize the variation of the complex momentum 

TTii'JirbUi multiplied by the factor rl{S/C), a complexifyed version of the planar force. 

To include the self-interaction into the linearized Euler-Lagrange equations, we simply add 

the complexifyed planar force of Eq. (10) multiplied by the same factor rl{S/C) into Eq.(41) 

9 (7*^7/ 29^ 

rl{S/C)-{l + mn-^ ^ ^(1 + 2^)-^i, (43) 

14 



where the dots represent derivative respect to the scaled time t = Qt. Using definition 
(11) for the gyroscopic coordinate ui into Eq. (43) yields the offensive force of Eq. (79) 
of Appendix C plus the required linearized gyroscopic version of the triple dots term. The 
offensive force is a nonhomogeneous term that will be dealt with in Appendix C. In this 
section we discard it and keep only the linear correction to the self force. 

Addding the linear term of the self-force to the rji, ^2 and 7)2 Euler-Lagrange equations 
yields a system of four linear delay equations, which can be solved in general by Laplace 
transform [15]. Here we shall concentrate on the separable normal mode solutions of this lin- 
ear system. For that we substitute = Aex.p{\T/6), rj^ = Sexp(Ar/6'), {2 = C'exp(Ar/^) 
and = Dexp{XT/6) into the linearized equations and we assume that |A| is an order-one 
or larger quantity, which is henceforth called the stiff limit (we shall see that |A| can take 
values of about tt or larger while 9 is about 10~^ in the atomic magnitude). The four linear 
equations in A, B, C and D have a nontrivial solution only if the determinant vanishes, 
a condition that we evaluate with a symbolic manipulations software in the large-A limit 
(stiff-hmit) 



In the special case of k — —1/2 Eq. (44) reduces to Eq. 33 of reference [2] without the 
radiative terms, that are ad- hoc in the dissipative Fokker setting of Ref.[2] (this amounts 
to setting = in the first line of Eq. (44) and k — —1/2 ui the other two lines). In this 
work we shall focus on Dirac's retarded-only electrodynamics of point charges (Eq. (44) 
with A; = 0) , with the following planar-normal-mode condition in the stiff-limit 



Comparing Eq. (45) to Eq. (78) of Appendix B we find that the quasi-degeneracy of the 
stiff dynamics exists only for A; = 0, (retarded-only interactions, Dirac's theory). A; = — 1 
(advanced-only interactions), and k — —1/2 (either the dissipative Fokker theory of Ref.[2] 
or the action-at-a-distance electrodynamics). 




(44) 




~ + •••)(iri-)exp(-2A,;j,) 



(45) 
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VI. STIFF TORUS AS A DEFORMED CIRCULAR ORBIT 



As discussed in the previous Section, in Appendix B and already explored in Ref.[2], there 
is a remarkable quasi-degeneracy of the perpendicular and the planar tangent dynamics for 
some special settings. Here we shall focus on Dirac's electrodynamics with retarded-only 
interactions. Both Eq. (78) and Eq. (45) with k — —1/2 in the large-A limit reduce to 

(^) exp(-2A) = 1, (46) 

For hydrogen {/jl/M) is a small factor of about (1/1824). For 9 of the order of the fine 
structure constant the small parameter ^ ~ 10~^^ multiplying the exponential function on 
the left-hand side of Eq. (46) determines that Re(A) = -cr ~ -ln{.^J^). For the first 13 
excited states of hydrogen this a is in the interval 14.0 < \a\ < 18.0. The imaginary part of 
A can be an arbitrarily large multiple of tt, such that the general solution to Eq. (46) is 

A = -C7 + Tiqi (47) 

where i = and q is an arbitrary integer. Notice that the real part of A is always 

negative, such that the tangent dynamics about the circular orbit is stable in the stiff-limit. 
The coefficient of the term of order 1/A^ is the main difference between Eqs. (78) and 
(45). The terms of order also differ but they are much less important for 6 in the 
atomic range. The exact roots of Eqs. (45) and (78) near the limiting root (47) are defined, 
respectively, by 

Xxy{0) = -a^y + Trqi + iei, (48) 
X;,{6) = -a^ + nqi + ie2, 

where ei(^) and €2(0) are real numbers of the order of 9 and depending on the orbit through 
9. The stiff torus is formed from an initial circular orbit as follows; (i) The offending force 
against the velocity dissipates energy and makes the electron loose radius on a slow timescale 
(the radiative instability discussed in Appendix C). (ii) When the electron deviates enough 
from the circular orbit, the stiff nonlinear terms compensate the offending terms (iii) Because 
of this same deviation from circularity, the stiff terms also balance the real parts of Eq. (48) , 
the a's . A purely harmonic solution to the nonlinear equations of tangent dynamics appears 
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because the nonlinear stiff terms introduce a correction into Eqs. (45) and (78), as discussed 
in Appendix D. The complete description of this new solution involves expanding the action 
about the circular orbit to the fourth order and the coupling to the center-of-mass recoil and 
shall be given elsewhere. Because of the large frequency of the stiff motion, this balancing 
is established at a very small radius, as estimated by Eq. (89) of Appendix D. The motion 
continues to be stiff because of the imaginary part of the X's of Eq. (48). We henceforth 
assume that this balance is achieved near the original circular motion. As discussed in 
Appendix C, the circular orbit recoils in a slow timescale to compensate the momentum 
loss to the radiated energy. During the recoil the particle coordinates are described by a 
composition of a translation mode and a stiff mode, both with slowly varying amplitudes, 
such that the dynamics of the balanced stiff torus is described by 

= Ak{T) + Uk{T) exp[(7rgi + ie>[)nt/ei (49) 
Zk = Re{B„{T) + RkiT)exp[{7iqi + ie'2)nt/9]}, 

where the amplitudes Uk and Rk must be near the size where the nonlinearity balances the 
negative real part of the linear modes. Assuming ~ ~ p, with p given by Eq. (89), 
the angular momentum of the stiff torus along the orbital plane calculated using (49) and 
disregarding fast oscillating small terms is 

k + tly = l^rlp^^h exp[2(e^ - e? + e)nt/e] (50) 
u 

= '^p^eM^{e'2-e1 + e)^lt/e], 

where on the second line of Eq. (50) we have used formula (1) for Q and formula (2) for rf,. 
Notice that we chose the same q on the two perpendicular oscillations in Eq. (49), such that 
the fast oscillations beat in the slow timescale. A typical nonlinear term of this dynamics is 
the angular momentum of spin defined by Eq. (49). As illustrated in Fig. 1, the balancing 
of the fast dynamics, Eq. (49), defines an angular-momentum vector (50) associated with 
the fast spinning. This angular momentum vector rotates at a (slow!) frequency that is 
generically of the order of the orbital frequency and is determined solely by the balancing 
of the fast dynamics. This slowly rotating gyroscope attached to the electron is the main 
qualitative feature left after we balance the stiff delay dynamics. According to Eq. (50) this 
slow frequency is 

+ (51) 
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Eq. (51) defines a frequency of the order of Q because 63 and are of order 9 by the 
definition of expansion (48). The spin angular momentum vector carries inertia, and its 
rotation frequency must influence the orbit, such that the slow trajectory should display 
some oscillation at the frequency zu of Eq. (51). We have assumed that the slow orbit 
is a circular orbit of frequency Q, therefore we must have that Q = zu \ (this is the 
simple physical consequence of solving for the fast dynamics first!). This physical resonance 
condition (Q = ro), by use of Eq. (51), yields 

- = (52) 

The calculation of and 63 necessitates and expansion of the action to fourth order and shall 
be given elsewhere. We expect nevertheless that the quasi-dcgcncracy is preserved, because 
the stiff limit depends basically on the time light takes to travel between the particles. Once 
the stiff torus is near the circular orbit, el and 62 should differ from ei and 62 by a correction 
which is essentially the same because of the quasi-degeneracy plus a correction of order 9 

ei - 62 ~ - + b9 (53) 

Condition (52) determines an orbital frequency proportional to a difference of two eigenval- 
ues, a Rydberg-Ritz-like formula 

O = 0(ei- 62)76^ = /x-(ei- 62). (54) 



where we have used Kepler's law of Eq. (1). The exact calculation of b demands expanding 
the action to fourth order and shall be given elsewhere. Since b must be of order one, we 
shall henceforth set 6 = —1 into Eq. (53) as a qualitative approximation, yielding 

61-62 + 9 = 0, (55) 

which has the solutions listed in Table 1. In Ref. [2] we had postulated heuristically a 
similar type of resonance for an clcctromagnctic-likc problem; one involving the real parts 
of the X's. The resonant orbits found in Ref. [2] were also in the atomic scale and had 
the same qualitative behavior found here, which seems to be a generic feature of these stiff 
electromagnetic resonances of the two-body problem. In the last paragraph of Appendix C 
we derive this same resonance condition (52) from the solvability of a multiscale asymptotic 
solution. The above-defined root-searching problem of Eq. (55) is well posed and for each 

18 



integer q it turns out that one can find a pair of the form (48) if one sacrifices 9 in Eqs. (78) 
and (45), i.e., 9 must be quantized! According to QED, circular Bohr orbits have maximal 
angular momenta and a radiative selection rule ( A/ = ±1) restricts the decay from level 
/c + 1 to level k only, i.e., circular orbits emit the first line of each spectroscopic series 
(Lyman, Balmer, Ritz-Paschen, Brackctt, etc.). the third column of Table 1. We have 
solved Eqs. (55) and Eqs. (78) and (45) with a Newton method in the complex A-plane. 
Every angular momentum 1/9 determined by Eq. (55) has the correct atomic magnitude. 
These numerically calculated angular momenta 9~^ are given in Table 1, along with the 
orbital frequency in atomic units (137^0)//^ = 137^^^(e2 — ei), and the QED first frequency 
of each spectroscopic series. 



h = 9-' 


1373^'(e2-ei) 


wqed 


Q 


185.99 


3.996x10-1 


3.750x10-1 


7 


307.63 


8.831x10-1 


6.944x10-2 


9 


475.08 


2.398x10-2 


2.430x10-2 


11 


577.99 


1.331x10-2 


1.125x10-2 


12 


694.77 


7.667x10-=^ 


6.111x10-3 


13 


826.22 


4.558x10-=^ 


3.685x10-3 


14 


973.12 


2.790x10-3 


2.406x10-3 


15 


1136.27 


1.752x10-3 


1.640x10-3 


16 


1316.44 


1.127x10-3 


1.173x10-3 


17 


1514.40 


7.403x10-3 


8.678x10-^ 


18 


1730.93 


4.958x10-3 


6.600x10-^ 


19 


1966.77 


3.379x10-^ 


5.136x10-^ 


20 


2222.70 


2.341x10-^ 


4.076x10-^ 
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Caption to Tabic 1: Numerically calculated angular momenta = 9~^ in units of e^/c, 
orbital frequencies in atomic units (137^)3 = 137^9'^{€2 — ei), circular lines of QED in atomic 
units Wqed = lip ~ (fe+Tp') ' values of the integer q of Eq. (48) and Eq. (49). 

Table 1 illustrates the fact that a resonance involving the tangent dynamics of a circular 
orbit predicts magnitudes in the atomic scale, as first discovered in Ref. [2]. In Ref. [2] 
we had to jump the integer q by twenty units for a complete quantitative and qualitative 
agreement with the Bohr atom. The qualitative agreement achieved by Table 1 is superior 
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in this way; only for the first two values of q things are not completely natural, and we had 
to jump some q's , after that q increases one by one in qualitative agreement with QED. 
It remains to be seen if using the correct b given by the expansion about the stiff torus 
will change things quantitatively. The calculation of b necessitates expanding the action to 
fourth order and shall de done elsewhere. 



VII. CONCLUSIONS AND DISCUSSION 



In the limit where the proton has an infinite mass, the concept of resonant dissipation 
looses meaning because the center-of-mass coordinate no longer plays a dynamical role. 
In this singular limit, there is a Lorentz frame where the proton rests at the origin at all 
times, and the field at the electron reduces to a simple Coulomb field in the ES. The two- 
body dynamics in the ES reduces then to the dynamical system of Ehezer's theorem; self- 
interaction plus a Coulomb field acting on the electron [5,6]. We repeat this correct dynamics 
because it is very unpopular [5-8]; With inclusion of self- interaction, it is impossible for the 
electron to "spiral into the proton". Neither bound states nor dives are possible, only 
scattering states exist. One accomplishment of the present work is to recognize that only 
the two-body problem can produce a physically sensible electromagnetic-like model. By 
expanding the dynamics about a quasi- planetary orbit (the circular orbit) , we saw the need 
for a novel ingredient of the electromagnetic two-body dynamics; the stiff spinning motions! 

Another qualitative dynamical picture is suggested by Eliezer's result [5, 6]; The dynam- 
ical phenomenon that the electron always turns away from the proton along unidimensional 
orbits suggests that colinear orbits are the natural attractors of the dissipative dynamics 
(the ground state of the hydrogen atom, with zero angular momentum!). Along such orbits, 
the heavy particle (the proton) moves in a non-Coulombian way and the self-interaction 
provides the repulsive mechanism that avoids the collision at the origin. This is again in 
agreement with the Schroedinger theory, where the ground state has a zero angular mo- 
mentum. Again, the infinite-mass case produces unphysical dynamics; the electron turns 
away from the proton but then it runs away [6]. It remains to be researched if the two-body 
version of Eliezer's problem has a physical orbit among its zero-angular momentum orbits, 
by taking proper care of the delay and its associated fast dynamics. 
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The angular momentum vector of fast spinning is rotating at the orbital frequency, such 
that the total angular momentum executes precession about an axis that is determined by 
the initial condition. Because the angular momentum of fast spinning is of the order of the 
orbital angular momentum, the orbital plane may also oscillate and nutate about an axis 
that is determined by the intial condition. The electromagnetic fields generated by such 
motion have this axis as a symmetry axis, and therefore this must be the axis of recoil, 
since the net radiated momentum calculated with the Poyting vector must be along the 
symmetry axis. The unbalanced momentum along the symmetry axis is the precise cause of 
the recoil. We expect that our planar circular orbit is a qualitative approximation to this 
complex slow motion. The correct multiscale expansion should not start from a circular 
orbit, but rather from a general slow motion of both particles plus a fast perturbation, and 
then balance the fast dynamics first. The linearized dynamics about these slowly moving 
positions is controlled by a 6 x 6 linear matrix containing the stiff delay modes, as the 
nonlinear terms couple the z and xy oscillations. This multiscale description should yield 
differential equations for the slow guiding orbit by a Preedholm alternative[25], and our 
preliminary findings indicate that these complex guiding orbits shall be found in the atomic 
region. This detailed research is yet to be done. 

In the dynamical process of resonant dissipation, the sharp line is emitted while the 
dynamics is locked to the neighborhood of the resonant orbit for a long timescale. This long 
timescale is to be compared to the time of spontaneous decay prescribed by QED for the 
circular hydrogen hues; about 10^ orbital turns or 10~^° seconds. This is indeed a very long 
timescale compared to the orbital period and should not be confused with the time for a stiff 
jump. The stiff modes of Eq. (48) have a frequency of about yrg/^ ~ 1000 times the orbital 
frequency. This fast frequency is about 10^° Hertz and resonates with the X-ray frequencies 
used in the Compton effect [18]. This also sets the timescale for a stiff jump of the dynamics; 
10~^^ seconds (10 attoseconds) . It is interesting to notice that this stiff frequency is exactly 
of the same magnitude of the zitterbewegung of Dirac's relativistic version of Schroedinger's 
equation[19]. In reference [19] it is also proposed that the spin should be associated with a 
dynamical motion. 

The dynamics starting from an asymptotic resonant orbit to another of a neighboring 
q should be described by a stiff jump, as expected generically from any stiff equation. In 
Ref [20], the much simpler Van der Pol oscillator is worked out in detail as an example of 
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an equation of Lienard type that exhibits stiff jumps. It is popular to use the concept of 
a quantum jump to describe the stiff passage from one quantum state to another, but the 
fact that classical electrodynamics prescribes exactly this qualitative phenomenon is news. 
This quasi-instantaneous fast dynamics could be calculated theoretically and compared to 
experiment. 

The two-body dynamics in the ES solves several conundrums of the hydrogen atom. Most 
of those conundrums were created by imagining that the equations of electrodynamics would 
accept nonstiff planetary-like orbits. We have seen that the stiff spinning is necessary if a 
balance of the fast dynamics is to be achived in the neighborhood of quasi-planetary orbits, 
as proved by the Lemma of resonant dissipation of Appendix C. The balanced stiff spinning 
is a novel and non-planetary feature that introduces a gyroscopic torque in the dynamics. 
The qualitative agreements with QED are listed in the following; (i) the resonant orbits are 
naturally quantized by integers and the radiated frequencies agree with the Bohr circular 
lines within a few percent average deviation, (ii) the angular momenta of the resonant orbits 
are naturally quantized with the correct Planck's constant, (iii) the stability analysis near 
the balanced stiff torus defines a linear dynamical system with delay, a dynamical system 
that needs an initial function as the initial condition, just like Schroedinger's equation. It 
remains to be seen if this linear operator produces a self-adjoint Freedholm alternative[25] , 
like in Schroedinger's equation, (iv) The emitted frequencies are given by a difference of two 
eigenvalues of this linear operator, like the Rydberg-Ritz combinatorial principle of quantum 
physics, (v) the spin angular momentum of the fast toroidal motion is of the order of the 
lowest orbital angular momentum. The spin that we estimate is still a bit too high for QED, 
where the electron has a spin angular momentum of \/3/i/2 , and perhaps the full equations 
can cure this[26, 27]. The calculation of the spin angular momentum for the balanced stiff 
dynamics involves a detailed consideration of all stiff terms and is beyond the estimates of 
the present work. 

We exhibited a new solution of two-body motion in Dirac's electrodynamics of point 
charges. The stiff dynamics appears naturally in the two-body dynamics because of the 
delay, and it has been so far overlooked. The balancing of the fast dynamics leads naturally 
to fast spinning motions; the multiscale analysis imposes resonances, and these turn out to be 
satisfied precisely in the atomic magnitude! The large body of qualitative and quantitative 
agreement with QED suggests that further extensive studies of this two-body dynamics [4] 
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could offer an explanation of QED in terms of a stiff dynamical system with third derivatives 
and delay. 
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IX. APPENDIX A: ELIEZER'S ELECTRODYNAMICS OF POINT CHARGES 

In Eliezer's generalized electrodynamics [4], the field produced by the point charge is 
supposed to be the retarded field plus a free field G 

F; = F;^ret + Gl- (56) 

The free field G should satisfy Maxwell's equations, be finite along the particle's worldline 
and vanish if the particle is at rest. The choice of Eliezer in [4] is 

Gl = k{F;^ret-Flavd). (57) 

where /c is a fundamental parameter of nature related to the covariant limit producing the 
point charge [4]. This generalized electromagnetic setting is henceforth called the Eliezer 
setting (ES). Analogously to Dirac's theory [1], Eliezer's self-interaction is given by the 
sourceless combination of half of the retarded Lienard-Wiechert self-potential minus half 
of the advanced Lienard-Wiechert self-potential, but multiplied by a renormalizing factor 
of (1 -I- 2k). This gives the following concise description of the ES; Charges interact with 
themselves via the semi-difference of Lienard-Wiechert self-potentials and with other charges 
via a linear combination of Lienard-Wiechert potentials. In Eliezer's theory the electron and 
the proton of a hydrogen atom have the following equations of motion [4] 

m,v,, - ^(1 + 2k)[v,, - \\v,\\\,,] = + (1 + /c)F;2,,et - kF^^Mvt]^^^^ (58) 

m2V2t.-\{l + 2k)[v2^,-\\v2\\\2^.\ = [i^;,,„ + (l + A;)F;i_,,,-A;F;i,„,>2., 

where Vifj^ stands for the quadri velocity of particle i, double bars stand for the Minkowski 
scalar product and the dot represents derivative respect to the proper time of each particle. 



23 



We are using units where the electron and the proton have charges —1 and 1 respectively and 
the speed of hght is c = 1. The ES has three important hmits; (i) For k — the ES reduces 
to Dirac's electrodynamics with retarded-only fields [1]. (ii) For k = —1/2 the ES reduces 
to the action- at-a-distance equations of motion derived from Fokker's Lagrangian [17] ( 
notice that the self- interaction terms disappear), (iii) For k ~ —1/2 the ES approximates 
the dissipative Fokker electromagnetic setting of Ref. [2] with a charge renormalization 
controlled by (1 -|- 2k). The ES is discussed in the excellent review of Ref. [4]. 



X. APPENDIX B: LINEAR STABILITY ANALYSIS ALONG THE z DIRECTION 

In this appendix we perform the linear stability analysis of the circular orbits for dis- 
placements perpendicular to the orbital plane, henceforth called the z-direction. We expand 
to second order the implicit light-cone condition and the action (7), in the same way of Sec- 
tion V. The linearized- 2; dynamics is uncoupled from the planar dynamics. The Cartesian 
coordinates of a transversely perturbed circular orbit are defined as 

Xk + iyk = ndk exp{iQt) , (59) 

Zk = fbZk, 

where k = 1 for the electron and k = 2 for the proton, Zk is the small transverse pertur- 
bation, di = bi and o?2 = — ^'2 are defined in Eq. (3) and fl is the orbital frequency defined 
above Eq. (3). We introduce again the delay function (p of the Zi and Z2 perturbations by 
expanding the light-cone time t2 about the constant lag rj, precisely by Eq. (15). In the 
following we calculate this homogeneous function ip of Zi and Z2 up to the quadratic order. 
The distance entering Eq. (9) is to be evaluated from the position of particle 1 at time 
ti, to the position of particle 2 at time t2 defined by Eq.(59). Using t2 defined by Eq. (9) 
and the orbit of particle 2 of Eq. (59), the implicit distance ru — \t2 — ti\ from particle 1 
at time to particle 2 at time ^2 is 

ru = rl{l + 4>f = rl[hl + hl + 26162 co&{^ + Oc)] + rl{Z^ - Z2cf . (60) 

Where we again expressed the </? of Eq. (15) in terms of the scaled function defined in 
Eq. (26). Notice that the Z variations decouple from the planar variations because there 
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is no mixed linear term of Z times a linear perturbation of the planar coordinate in Eq. 
(60); terms in Z only appear squared. Expanding Eq. (60) up to the second order and 
rearranging yields 



02 + 250 = (Zi - Z2e)' (61) 

a quadratic equation for with the regular solution correct to second order given by 

= ^(^1 - Z^cf- (62) 

The coordinate Z2 appears evaluated at the advanced/retarded time in Eq. (62), and to 
obtain the action up to quadratic terms it is sufficient to keep the first term Z2{t\ + + (/?) 
— Zi{t\ + c^) = Z2c- Using the perturbed orbit defined by Eq. (59) to calculate the 
numerator of the VA interaction of Eq. (8) yields 

(1-Vi-V2c) = 1 + ^2 C0S(^ + C(^)6i62 - ^^^1^2c ~ (63) 

C - B^{S - 1)0 - Q^Z^Z^c, 

and the denominator of the VA interaction of Eq. (8) is 

ri2(l + ni2c • V2c/c) = rfe[l + + Qchxh^ sin(^c + (^J + Qc^Z^ - Z2c)Z2c]. (64) 

Notice that the quadratic term Z2cZ2c on the right-hand side of Eq. (64) can be dropped 
because it represents an exact Gauge that does not affect the Euler-Lagrange equations of 
motion. We also expand the argument of the sign function of the right-hand side of Eq. (64) 
until the linear term in <^^, such that the quadratic approximation to Eq. (64) is 

ri2(l + ni2c • V2c/c) ^ n\S + C0 + QcZ^Z2^, (65) 

where the equivalence sign henceforth means equivalent up to a Gauge term of second 
order. Even if a quadratic Gauge term appears in the denominator, in an expansion up to 
quadratic order it would still produce a Gauge and therefore it can be dropped directly from 
the denominator. In this way, the expansion up to second order of the VA interaction of Eq. 
(7) is simply 

e> ~ (-^){(i + fc)[i - Q''cs''z^Z2- - ^{Zi - Z2-f + ec''sz^Z2-] (66) 

-A;[l - e''CS^ZiZ2+ - ^{Zi - Z2+f - eC''SZiZ2+]}. 

25 



Last, we need the kinetic energy along the 2;-perturbed circular orbit, which we express in 
terms of Z\ of definition (59) as 

^1 = -mi0;^= -^^1 - ^\Q^Zl (67) 

where the dot means derivative with respect to the scaled time r, 7f ^ = ^Jl — v\ , and we 
have used Qr^ = d. The expansion of Eq. (67) up to second order is 

= (^){^ + |Z? + ■..}. (68) 

where ei = m\r\pi^Q^ — r\m\'^^ . The Euler- Lagrange equation of motion for particle 1 of 
the isolated two-body problem is determined from the quadratic Lagrangian 

Li = Ti + e. (69) 

Wc shall henceforth disregard the small corrections to the quadratic coefficients and use 
C = 1, 5* = 1 , which is the stiff-limit. The equation of motion for particle 1 is 

eiZi = -[Zi - (1 + A;)Z2_ + kZ^^] + Q\kZ2^ + (1 + fc)Z2_) + Q^kZ^^ - (1 + k)Z-,_\. (70) 

Notice that the term on the left-hand side of Eq. (70) can be written as 

6iZi = rlm^^.Vi'Z^ = r,^^, (71) 

which is proportional to the force along the z-direction multiplied by the factor r^. According 
to the equation of motion of the ES (see appendix A), we must add the following self- 
interaction term to the right-hand side of Eq. (70) 

r,2F,„rf = ^(l + 2A;)^^Zi, (72) 

where the triple dot means three derivatives with respect to the scaled time and we have 
used Eq. (10). The full hnearized equation of motion for Z\ is 

eiZi = ^(1 + 2A;)^'Zi - [^i + kZ.,^ - (1 + k)Z.,_\ (73) 
+^[A;Z2+ + (1 + k)Z2-\ + Q\kZ-,^ - (1 + k)Z2-\ 

The linearized equation for Z^ is completely analogous and is obtained by interchanging 
Zi by Z2 and ei by e-i in Eq. (73). (Comparing Eq. (73) to Eq. (30) of Ref. [2] we find that 
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Eqs. (29) and (30) of Ref.[2] are both missing a 9^ factor in front of the Zi term, which 
is just a typo in Ref. [2] because from Eq. (31) on in Ref [2] the self-interaction force 
was included correctly). The general solution of a linear delay equation can be obtained by 

Laplace transform [15] and is a linear combination of the following normal mode solutions. A 
normal mode solution is obtained by substituting Zi = Aexp(pr) and Z2 = -Bexp(pr) into 
the two linearized equations, and requires the vanishing of the following 2x2 determinant 

l + eip2-|(i + 2A;)eV G{e,p) 

G{e,p) l + ea/- |(l + 2A;)^V 

where G(^,p) = [1- {l + 2k)p9 -9'^p'^]cosh{p9) -[{l + 2k) -9p- {l + 2k)9'^p'^]smh{p9). The 
stiff limit obtained when p9 is large, such that the hyperbolic functions of the G{9, p) acquire 
a large magnitude [16]. In the following we use the Coulombian limit values, hi — 1712/ M 
and 62 = rrii/M , (see Appendix B of Ref. [2]) to evaluate determinant (74), such that 



detZ 



(74) 



ei = — , (75) 

1712 

M 

£2 = ■ 

nil 

For small 9, the second-order and higher even-order corrections to Eq. (75) are very small. 
Defining p = X/9 and using Eq. (75) we obtain 

(detZ) = i-l(i + 2k)9^X+--^9^X'' (76) 



-{[1 - i + ^^^] cosh(A) - [U (1 + 2k){l - h] sinh(A)}^ 



M X' X ' V / . V . .V ^2 

where we have dropped small 0{9^) terms. Equation (76) exhibits the generic feature 
of the stiff-limit; that the hyperbolic functions always appear multiplied by the very small 
coefficient ii9* /M . The stiff-mode condition defined by Eq. (76) (det Z = 0) with k — —1/2 
is equation (33) of Ref. [2], i.e., 

^ - l^'^ + - W^^^ - - ^sinh(2A)]^ = 0. (77) 

Notice that in Ref [2] there is a typo in passing from Eq. (33) to Eq. (34); Eq. (34) is 
missing a bracket that should start after the ^ factor and close at the end of Eq. (34). 
The special case of = is Dirac's theory with retarded-only fields 

' - 1"'' ^ m'''' - w^^^^M-2mi + ^ - ^ - J, , J,) = 0, (78) 
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where the appearance of the negative exponential only is related to the retardation-only, in- 
stead of the hyperbolic functions related to advanced and retarded interactions. Comparing 
Eq. (76) to Eq. (44) we can see that the phenomenon of quasi-degeneracy exists only for 
k = (Dirac's theory) and for k = —1/2 (the action-at-a-distance electrodynamics and the 
dissipative Fokker theory of Ref. [2]). 



XI. APPENDIX C: LEMMA OF RESONANT DISSIPATION 



Substituting the circular orbit into the equations of motion, Eqs. (58), one finds at 
leading order for the electronic motion an offending term described by a force opposite to 
the electronic velocity 

rlF, = rf-x, = = -le\ (79) 

where we have included the scaling factor as of Eq. (72). Along the trajectory of the 
proton, the delayed Lienard-Wiechert interaction with the electron is the main offending 
force against the velocity, instead of the much smaller protonic self-interaction. Using the 
Page series in the same way of Ref. [13], we find a force against the protonic velocity of the 
same magnitude of Eq.(79), as illustrated in Fig. 2. These offending forces destabilize the 
unperturbed circular motion with a small torque and a slow dissipation. In quantum elec- 
trodynamics (QED) the circular Bohr orbits[3] of hydrogen correspond to excited states that 
decay to the ground state in a life-time of about 10^ turns. Using Eq. (79) to estimate the 
energy dissipation along circular orbits of the atomic magnitude, one finds a net dissipation 
of about 4 electron-volts after the life-time of 10^ turns, an order-one fraction of the biding 
energy (13.6 electron-volts). This radiative instability in the slow timescale is an incomplete 
picture of the dynamics because it disregards the fast stiff motion that is present because 
of the delay. Let us postulate heuristically that along some special circular orbits the fast 
dynamics encircles the circular orbit as illustrated in Fig. 1. If the linearization about the 
circular orbit, Eq. (49), is averaged over a timescale of some turns, the stiff dynamics goes 
away and the resulting variable is a slow moving drift (representing the recoil of the state 
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of resonant dissipation). For that slow motion we have 

liir) = lMr)^l,_{r), (80) 
Ur) = e2+(r)=e2-(r), 

since we are averaging the argument over a time of many turns, a time much longer than 
the delay 6. We are ready to prove that the dynamics of the state of resonant dissipation 
must necessarily involve a nonlinear term of the expansion about a circular orbit. 

Lemma: If the averaged dynamics is determined only by the linear terms plus the offend- 
ing forcing, the state of resonant dissipation is impossible. To show this we average the full 
equation of motion of the variable, which is 

(—{ui + ibi))+h^-ihe'' + NL = Q, (81) 

where the term inside brackets is the average of the linearized correction to the self- 
interaction force, li^ stands for the average of Eq. (41) and NL stands for the average 
of the nonlinear terms. The time averages of the term inside brackets and of the first and 
second derivatives in Eq. (41) are zero, such that discarding the nonlinear term of Eq. (81) 
yields the non-homogeneous linear equation 

-(TVn + 0'G)C, - (Ml + C/n)77i (82) 

-iRl + R*_)l,-{P+ + P.)v, 

Performing the same average for the other three coordinates, we obtain a system of four 
non-homogeneous linear equations. The matrix of this linear system has a zero determinant 
that is directly related to the Galilean translation mode of the circular orbit, such that it 
defines a singular linear operator. Moreover, the forcing term turns out to be out of the 
image of the singular linear operator, such that the linear system has no solution at all! This 
completes the proof of our lemma. 

Last, we mention a derivation of Eq. (52) from the equations of motion in the form (81); 
This can be accomplished by placing the linearization about the stiff torus on the left-hand 
side (henceforth called the Kernel) and the forcing and nonlinear terms on the right-hand 
side. In this way the right-hand side has to be orthogonal to the left-eigenvectors of the 
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Kernel, which is the Preedholm alternative theorem discussed for example in Refs. [24, 25] . 
Multiplying the Kernel by a left-eigenvector and integrating on the fast timescale has exactly 
the form of the conservation law for angular momentum on the slow time scale. This shall 
be studied elsewhere. 

XII. APPENDIX D: HARMONIC SOLUTION AND SPINORIAL UNFOLDING 

The lemma of resonant dissipation of Appendix C suggests the participation of the nonlin- 
ear terms to achieve the state of resonant dissipation of Fig. 1. We shall see in the following 
that because of the fast frequency in the stiff terms, even at small amplitudes the nonlinear 
stiff terms give a large contribution to the self-force. Our perturbation scheme started from 
a guiding circular orbit that is not a solution to the equations of motion, because of the 
offending force of Eq. (79). In the following we show how to balance this offending force 
at a small distance away from the circular orbit. The full Lorentz-Dirac self-force is given 
in page 116 of Ref [22], expressed in terms of the Cartesian velocity of the particle, vi, its 
acceleration ai and its second acceleration ai respectively. This full self-force multiplied by 
the convenient scahng factor of Eq. (72) is 

rlFi = plrliki + 37'(vi • ai)ai + 7^[(vi • ai) + 37^(vi • a^Vi}. (83) 

Assuming that the particle coordinates are described by a slow term plus a fast term, Eq. 
(49), we use Eq. (83) as follows; We estimate each coefficient of the vi, ai and ai terms by 
substituting the fast component of the trajectory, Eq. (49), and taking the time average. For 
the estimate below, the vectors Vi, ai and ai are replaced by the slow quantities evaluated 
along the unperturbed orbit. For example the coefficient of the velocity term in Eq. (83) 
{rl times the term inside the square brackets on the right-hand side of Eq. (83) ) is 

rhi[(vi • ai) + 3^l{v^ ■ ai)^] = 7?|A| V(37?|A| V' - 1) (84) 

Notice that for a large enough p there is a critical value where this coefficient changes sign, 
or 

37?|A|V > 1, (85) 

We show in the following that the value of must be very near this critical value. The 
self-forces for the electron and proton, both measured along the direction of the electronic 
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velocity, are 



rlF, = ^H^&i + 7?|A|V?(37?|A|V? - (86) 



3' 

where hi and 62 are given by Eq. (3). In the state of resonant dissipation, the relative 



rlF, = l[e% - jl\X\y,{3jl\X\'pl - l)9h], (87) 



separation of the particles executes a circular motion, such that the tangential force must 
vanish 

hFi - = 0, (88) 

with bi = 1712/ {rrii + 1712) and 62 = rni/{mi + 7712). Because 62 ^ 1 we can disregard the 
protonic contribution such that Eq. (88) yields 



This distance is a few percent of the orbital radius, and it is some 500 classical electronic 
radia for 9 in the atomic scale. The spin angular momentum of the electron calculated by 
Eq. (50) is \l\ ~ l/(67rg^^), which is of the order of the orbital angular momentum. 

The nonlinear correction to Eq. (78) can be estimated by taking into account the non- 
linear terms of the Lorentz-Dirac self- force. This correction to Eq. (78), estimated along 
a harmonic solution for the averaged equations of tangent dynamics, in a way explained 
above, yields 



= 1 - 2p-'\\f - je^x, + ^e-Al + .... 

Because of the extra term on the right-hand side of Eq. (90), a harmonic solution to the 
nonlinear equations of tangent dynamics exists if 

p~ , . (91) 

The fact that the two necessary and independent estimates (89) and (91) agree suggests 
that the stiff torus of Fig. 1 is a solution of Dirac's equations of motion[l] for the two- 
body problem! In summary, our perturbative scheme shows that the stiff torus of Fig. 2 

31 



is a solution to the equations of motion as follows; (i) We expand about a circular orbit 
and calculate the radius of the stiff torus by the condition that a harmonic solution to the 
nonlinear equations of tangent dynamics exists (Eq. 91). (ii) This balanced stiff spinning 
about the circular orbit cancels the offending force of Eq. (79), predicting a radius by (89) 
that is in agreement with the radius of (91). (iii) Resonance condition (52) is a necessary 
condition for the existence of an asymptotic expansion of the dynamics about the stiff torus, 
i.e., the condition for a smooth recoil of the atom as a bound state. The detailed unfolding of 
the equations of tangent dynamics is beyond the estimates of this Appendix and necessitates 
the full nonlinear equations. 

XIII. CAPTIONS 

Fig. 1: The guiding-center circular orbit of each particle is illustrated by dashed hues. 
Particle trajectories are the stiff tori gyrating about the guiding-center circular orbit of each 
particle (dark solid lines). Illustrative purposes only, the radia are not on scale. Arbitrary 
units. 

Fig. 2: The unperturbed circular orbit with the particles in diametral opposition at the 
same time in the inertial frame. Indicated is also the advanced position of particle 1 and the 
angle travelled during the hght-cone time. The drawing is not on scale; The circular orbit 
of the proton has an exaggerated radius for illustrative purposes. Arbitrary units. 
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